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Effects of Dynamical Quarks in UKQCD Simulations 
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Kingdom 

Recent results from the UKQCD Collaboration's dynamical simulations are presented. The main feature of 
these ensembles is that they have a fixed lattice spacing and volume, but varying sea quark mass from infinite 
(corresponding to the quenched simulation) down to roughly that of the strange quark mass. The main aim of 
this work is to uncover dynamical quark effects from these "matched" ensembles. We obtain some evidence of 
dynamical quark efi'ects in the static quark potential with less effects in the hadronic spectrum. 



1. Introduction 

The UKQCD Collaboration has embarked 
upon a study of unquenching effects in QCD lat- 
tice simulations particularly in the static quark 
potential, light hadron spectrum and topologi- 
cal sector. The philosophy that UKQCD uses in 
this study is to simulate at a variety of sea quark 
masses, m**^", but with a fixed the lattice spacing, 
a. This is in order to maintain, as far as possi- 
ble, constant lattice systematic effects due to the 
finite-ness of a and the volume. 

It is well known that the lattice cut-off, a, is 
a function of both the gauge coupling, and the 
dynamical quark mass, m'"^" (see e.g. |^]). For 
this reason, the philosophy we have chosen is to 
simulate at points in the {l3,m^^°-) plane which 
have constant a. We term this the "matched" 
trajectory. The hope is that this procedure will 
disentangle lattice systematic effects (due to the 
finite-ness of a and volume) with unquenching ef- 
fects. Any variation of a physical quantity along 
this matched trajectory can be attributed to un- 
quenching effects rather than, e.g. 0{a^) effects. 

This paper presents an overview of this work. 
See for a full description. 



term was non-perturbatively determined by the 
Alpha Collaboration Q. 

A range of {(3, k^'^'^) values were chosen in or- 
der to maintain a constant value of the lattice 
spacing. This used the technology of ji). The pa- 
rameter values for the simulations are displayed 
in Table |^. The lattice volume used was 16'^ x 32 
throughout. 

Table Q lists also the lattice spacing obtained 
from the Sommer scale, The first error 

listed is statistical, and the second (shown as tl) 
is the systematic error from variations in the fit 
used for the static quark potential. The central 
values quoted were obtained using all potential 
data satisfying V2 < r < 8. It can be seen 
that the last four simulations listed (i.e. those at 
= 0.1350,0.1345,0.1340 and 0) are matched 



to within errors, with the 



= simulation 



corresponding to a quenched run. 

The simulations at k"'^^ = 0.1355 and 0.13565 
were performed in order to study lighter sea quark 
masses and do not lie on the above matched tra- 
jectory. 

Full details of autocorrelation times and other 
algorithmic issues can be found in 1^. 



2. Simulation Details 

The standard Wilson gauge action was used 
together with the Clover C'(a)-improved fermion 
action. The coefficient, cswi used in the Clover 



3. The static potential 

The static quark potential was determined for 
all our ensembles using the method originally 
proposed in and using a variational basis of 
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Table 1 

Lattice parameters together with measurements of a and quark masses. 
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Figure 1. The static QCD potential expressed in 
units of rg. The dashed curve is a string model 
described in the text. 



Overall there is no discernable difference be- 
tween the ensembles at distance r « rp. Fur- 
thermore, the data follow the string model very 
well. However, at shorter distances, r < 0.5ro, 
there is a deviation from the string model, and 
a variation amongst the ensembles. This seems 
to be systematic in the sea quark mass - the de- 
viation from the string model increases as to**^" 
decreases. In fact, careful correlated fits of the 
potential show that the parameter e in Eq. (|lj) in- 
creases by 181io% ill going from the quenched to 
^soa _ .13500 data. We note that the ensembles 
in Figure 2 are matched, and therefore we can ex- 
clude lattice systematics from the effects that we 
have observed. (Similar findings in the case of two 
flavours of Wilson fermions have been reported 
by the SESAM-TxL collaboration g where an 
increase of 16 — 33% was found.) Finally we note 
that there is no evidence of string breaking in 
our static quark potential. However the distance 
scales covered r ^ 1.3 fm is not large. 



Figure 1 plots the static potential in units of 
To . The zero of the potential has been set at r = 
tq. The data are well described by the universal 
bosonic string model potential [^ which predicts 

[V{r)-V{ro)]ro = (1.65 - e) (^^ - 1^ 

Of course, the fact that the scaled potential mea- 
surements all have the same value and slope at 
r = ro simply reflects the definition of rg. Figure 
2 shows the deviations from this model potential 
for the matched ensembles. Here e = 7r/12 [|j. 



4. Hadronic Spectrum 

In this section one of our main aims will be to 
uncover unquenching effects in the light hadron 
spectrum. Because we have a matched data set, 
any variation amongst our ensembles can be at- 
tributed to unquenching effects. However, the 
task of identifying variations is likely to be hard 
for those quantities which are primarily sensitive 
to physics at the same scale as that used to define 
the matching trajectory in the (/3, k^'^'^) parame- 
ter space (ro in this case). This is expected to be 
the case for the hadron spectrum considered here 
where the quark masses are still relatively heavy. 



3 



• 5.93 Quenched 

5.29, C.I .92, k=. 13400 
□ 5.26, C.I .95, k.. 13450 
- 5.20, c.2.02, k=. 13500 



Figure 2. The difference between the static QCD 
potential expressed in physical units and the pre- 
diction of the string model described in the text. 
For clarity, only data from the matched ensembles 
are shown. 



Two-point hadronic correlation functions were 
produced for each of the ensembles using interpo- 
lating operators for pseudoscalar, vector, nucleon 
and delta channels as described in Mesonic 
correlators were constructed using both degener- 
ate and non-degenerate valence quarks, whereas 
only degenerate valence quarks were used for the 
baryonic correlators. 

In the following, we review the main fitting 
procedures which were used to obtain the light 
hadron spectrum results. Further details of the 
fitting procedure can be found in 

We used the fuzzing procedures of to gen- 
erate correlators of the type LL, FL and FF where 
F denotes fuzzed, and L local operators. 

Correlated fits were used throughout the fit- 
ting analysis of the correlation functions with the 
eigenvalue smoothing technique of of [ p^ em- 
ployed. A factorising fit was performed which 
combined the three fuzzed correlators LL, FL and 
FF together to extract a better estimate of the 
ground state parameters. 

Full details of the fitting procedure can be 
found in B. 



4.1. The J parameter 

In Figures 3 and 4 the vector meson masses and 
hyperfine splittings are plotted against the corre- 
sponding pseudoscalar masses for all the datasets. 
It is difficult to identify an unquenching signal 
from these plots - the data seem to overlay each 
other. Note that in it was reported that there 
was a tendency for the vector mass to increase 
as the sea quark mass decreases (for fixed pseu- 
doscalar mass). The observations for the present 
matched dataset imply that this may have been 
due to either an 0{a) effect (since the dataset in 
10] was not fully improved at this level) or a finite 
volume effect. The conclusion therefore is that it 
is important to run at a fixed a in order to disen- 
tangle unquenching effects from lattice artifacts 
or finite volume effects. 
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Figure 3. Vector mass plotted against pseu- 
doscalar mass squared in units of tq, together 
with the experimental data points. 



A possible explanation for the lack of convinc- 
ing unquenching effects in our meson spectrum 
is the following. Our matched ensembles are de- 
fined to have a common rg value, so any physical 
quantity that is sensitive to this distance scale 
(and the static quark potential itself) will also, 
by definition, be matched. Our mesons, because 
they are composed of relatively heavy quarks, are 
examples of such quantities, and this is a possi- 
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ble reason why there is no significant evidence of 
unquenching effects in the meson spectrum. 

A further point regarding hyperfine splitting in 
Figure 4 is that the lattice data for the matched 
ensembles tends to flatten as the sea quark mass 
decreases. (The quenched data has a distinctly 
negative slope, whereas the k^'^^ = 0.1350 data 
is flat.) Thus the lattice data is tending to- 
wards the same behaviour as the experimental 
data which lies on a line with positive slope (inde- 
pendent of the value used for tq). This behaviour 
is apparently spoiled by the unmatched run with 
^soa _ 0.1355 (see Figure 4) which has a clear 
negative slope. However, the k^'^'^ = 0.1355 data 
does not satisfy the finite volume bound of |^], 
and therefore we can attribute the trend in this 
data to a lattice artefact. 

We now study the J-parameter defined as pSl 



J = M 



dMv 



(2) 



K,K* 



along what we will term the 'unitary' trajectory, 
i.e. along k*'^^ = k™'. In Table ||, the results 
from both methods are given. These values of 
J are around 25% lower than the experimental 
value Jexpt = 0.48(2). 



/? 


^sca 


J 


First Approach 


5.2000 


0.1355 


0.321^ 


5.2000 


0.1350 


0.3931^" 


5.2600 


0.1345 


0.365i« 


5.2900 


0.1340 


0.349i^ 


5.9300 


0.0000 


0.3761^2 


Second Approach 









Third Approach 

0.43;^, 

Table 2 

J values from the various approaches as described 
in the text. 
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Figure 4. The hyperfine splitting. 



In the context of dynamical fermion simula- 
tions, this parameter can be calculated in two 
ways. The first is to define a partially quenched J 
for each value of the sea quark mass. In this case, 
the derivative in (H) is with respect to variations 
in the valence quark mass (with the sea quark 
mass fixed). The second approach is to define J 



Finally we note that the physical value of J 
(i.e. that which most closely follows the proce- 
dure used to determine the experimental value of 
Jexpt = 0.48(2)), should be obtained from extrap- 
olating the results from the first approach to the 
physical sea quark masses. We call this the third 
approach. In order to perform this extrapolation, 
we extrapolate the three matched dynamical J 
values obtained from the first approach linearly 
in {M'^f'^'Y to (M™'*"^'')2 = 0. M]if is the 
pseudoscalar meson mass at the unitary point. 
The value for J from the third approach is pre- 
sented in Table |^ and we note that it is approach- 
ing the experimental value for J . 

The results from all three approaches are plot- 
ted in Figure 5, together with the experimental 
result. There is some promising evidence that the 
lattice estimate of J increases towards the experi- 
mental point as the sea quark mass decreases (see 
the J value from approaches 1 and 3). 

Recently there has been a proposed ansatz for 
the functional form of My as a function of Mpg 
jl^ . Although this ansatz would be interesting 
to pursue, all our data have Mps/My 0.6, and 
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Experiment 



Approach 3 : 



Approach 2 



(Quenched) 



O Approach 1 



1131 0.132 0.133 0.134 0.135 0.136 0.1.37 0.138 0.L39 0.14 



Figiire 5. J versus k"^'^ using the approaches 
as described in the text. Note that the quenched 
data points have been plotted at k^°^ = 0.132 for 
convenience. Approaches 2 & 3 are obtained after 
a chiral extrapolation and are shown as banded 
regions. The experimental value J = 0.48(2) is 
also shown. 



for this region, the ansatz of is linear to high 
precision. Therefore we have chosen to interpo- 
late our data with a simple linear function and 
await more chiral data before using the ansatz of 

El- 

4.2. Lattice spacing 

This subsection presents a determination of 
a from the meson spectrum which complements 
that from tq. 

A common method of determining a from the 
meson spectrum uses the p mass. However, this 
requires the chiral extrapolation of the vector me- 
son mass down to (almost) the chiral limit which, 
as was discussed in the previous subsection, may 
be problematic. An alternative method of ex- 
tracting the lattice spacing using the vector me- 
son mass at the simulated data points (i.e. with- 
out any chiral extrapolation) was given in [p^ . 
Using this method, we obtain the lattice spacing 
values as shown in Table [| labelled aj . Note that 
these are in general 10-15% larger than the val- 
ues from Section || where the lattice spacing was 
determined from tq. A possible explanation for 
this discrepancy is that the potential and mesonic 



spectrum are contaminated with different 0{a?) 
errors (or that the value tq = 0.49 fm is 10-15% 
too small!). 

In order to investigate unquenching effects in 
the meson spectrum, we define the quantity 

<5,,(/3,m-) = 1^^4§^' 

where is the lattice spacing determined from 
the physical quantity i = {Mp, Mk, J-k ■ ■ ■}• Note 
that when 5i,j = 0, the lattice prediction of Af^ 
with scale taken from Mj agrees with experiment. 
Thus (5 is a good parameter to study unquench- 
ing effects. We expect that Sij{P,m'''^°-) = O(a^) 
since we are working with a non-pcrturbatively 
improved clover action. 



In Figure 6, 



is plotted against 



(aMp"*'*'^-*')"^ for the matched datasets. In this 
plot we have fixed j = vq and the various physical 
quantities i are ^/a (the string tension) and the 
hadronic mass pairs {Mk*,Mk) & {Mp^Mt^). 
The method used to determine the scale from 
these mass pairs is that of . It is worth noting 
that the experimental point on this same plot 
would occur at an x— co-ordinate (depressingly) 



of (aAfp 



unitaryx— 2 



200. 



O K/K< 
A p/7I 



l/{aM (nnitary))" 



Figure 6. 5i as a function of l/(aMp™*"'^^)^ 
for i = y/a and the mass pairs {Mk*, Mk) & 
{Mp, Mtt). Si is defined in eq.(||) with j = rp. 



Figure 6 does not show signs of unquenching 



6 



for quantities involving the hadronic spectrum 
i.e. the mass pairs {Mk*,Mk) & (Mp,M7r). 
However, there is evidence of unquenching effects 
when comparing the scale from rg with that from 
^/a since the quenched value of is distinct 
from the dynamical values. 

One may wonder if the the S values may have 
been distorted by not choosing the simulation pa- 
rameters {f3, m*'^") exactly on the matched trajec- 
tory. In order to obtain a rough estimate of the 
effect of mis-matched value of /3, we use the renor- 
malisation group inspired ansatz for a.^ 16|r^: 



a.(5o) = A- Vpt(5o) x [l + ^./pt(<7o)' 



(4) 



where fprig^) is the usual asymptotic scaling 
function obtained from integrating the /3-function 
of QCD and Xi is the coefficient of the ©(a") 
lattice systematic. The functional form for a(.gg) 
was originally applied for the quenched theory, 
but let us assume that it can also be applied in 
the unquenched case. Using Eq.(|), we sec that 
a mismatch in /? of A/3 would lead to a relative 
error in S of 



5{P + A/3) - 

m 



-3A/3. 



(5) 



This shows that even an error in f3 of as much as 
A/3 w 0.01 introduces a relative error in S{(3) of 
only 3%, ruling out the possibility that a signifi- 
cant distortion in 6 could have occurred due to a 
mismatching in /3. 

4.3. Edinburgh Plot 

In Table |l| the ratios Alps /My are displayed 
for the unitary case k^°^— k™^ As can be seen 
the simulated data is a long way from the experi- 
mental value M-jr/Mp = 0.18. Figure 7 shows the 
'Edinburgh plot' {Mn/Mv v.s. Mps/My) for all 
the data sets. There is no significant variation 
within the dynamical data as the sea quark mass 
is changed, but the dynamical data does tend to 
lie above the (matched) quenched data. This lat- 
ter feature may be indicative of finite volume ef- 
fects since these are expected to be larger in full 
QCD compared to the quenched case [Esl. 
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P = 5.20k=0.1355 
5.20k = 0.1350 
5.26 K= 0.1345 
5.29 K= 0.1340 
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_J 
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Figure 7. The Edinburgh plot for all the data 
sets. All degenerate k™' correlators have been in- 
cluded. The phenomenological curve (from fig] ]) 
has been included as a guide to the eye. 



5. Chiral Extrapolations 

In 1^] we used 3 approaches to perform the 
chiral extrapolations of hadron masses: "Pseudo- 
Quenched" ; "Unitary Trajectory" ; and a "Com- 
bined Chiral Fit" . In this publication we will refer 
only to the last approach. We take 

= y4(K'^°^) + S(k'''=^)Mps(k"°^; 
= Ao + AiMps(k''^«:''^)"^ 
+ \Bo -f BiMpsiK"""] K""")-^] Mpsin^"''; 



using the nomenclature M = aM. The first argu- 
ment of M{k^°^; k™') refers to the sea quark and 
the second to the valence quark. 

The results of these extrapolations are shown 
in Table ^. We stress that this functional form for 
the extrapolation is not motivated by theory, but 
is used as a numerical analysis technique in order 
to test for evidence of unquenching effects. As can 
be seen from Table ||, the parameters Ai and Bi 
are compatible with zero (to 2a) and therefore we 
conclude that there is no evidence of unquenching 
effects (i.e. there is no statistical variation in the 
fit parameters A and B with sea quark mass k^°'^). 
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hadron 



Bo 



Vector 

Nucleon 

Delta 



.663111 
.841^ 



-0.004tf, 
0.0061^ 
-0.002t? 



0.91+1 



.0151? 
-.0011} 
.0211 



Table 3 

Fit parameters from the Chiral Extrapolations 



6. Conclusions 

This paper attempts to uncover unqucnching 
effects in the dynamical lattice QCD simulations 
at a fixed (matched) lattice spacing (and vol- 
ume) and various dynamical quark masses. This 
approach allows a more controlled study of un- 
qucnching effects without the possible entangle- 
ment of lattice and unquenching systematics. 

We see some sign of unquenching effects in 
the static quark potential at short distance, but 
no significant sign of unquenching effects in the 
meson spectrum. This is presumably since our 
dynamical quarks are relatively massive, and so 
the meson spectrum is dominated by the static 
quark potential. This potential is, by definition, 
matched amongst our ensembles at the hadronic 
length scale tq, and so any variation of the me- 
son spectrum within our matched ensemble must 
surely be a "higher" order unquenching effect 
which is beyond our present statistics. 

It is interesting to note that other work (using 
the same ensembles) has shown interesting un- 
quenching effects in the glueball and topological 
sector 01. 
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